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ABSTRACT 

A  nonlinear  control  system  sometimes  oscillate  at  a 
frequency  that  is  an  integral  submultiple  of  the  driving 
frequency  in  response  to  a  sinusoidal  input.   Such  a  re- 
sponse is  undesirable.   An  analog  computer  study  is  under- 
taken to  investigate  the  occurrence  of  subharmonic  oscillation 
in  second  order  systems  with  nonlinearities  characterized  by 
saturation  and  hysteresis  effects. 
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I.   INTRODUCTION 

A.    NONLINEAR  SYSTEMS 

A  system  is  said  to  be  linear  if  its  performance  obeys 
the  "Principle  of  Superposition".   The  "Principle  of  Super- 
position" states  that  if  C, (t)  and  C?(t)  are  the  responses 
of  a  system  when  it  is  separately  subjected  to  two  inputs 
R-.  (t)  and  R~(t),  respectively,  then  for  all  a  and  b,  the 
response  of  the  system  to  the  input  aR. (t) +bR- (t)  is 
aC,  (t)  +bC;?  (t)  and  this  must  hold  for  all  inputs.   In  other 
words,  once  the  response  to  one  type  of  signal  is  known, 
the  response  to  any  other  signal  can  be  deduced  directly. 

On  the  other  hand,  any  system  is  said  to  be  nonlinear 
if  the  "Principle  of  Superposition"  does  not  hold.   The 
term  nonlinear  systems  applies  to  varied  types  of  systems 
which  have  practically  nothing  in  common  with  one  another. 
As  a  consequence,  there  is  no  way  of  generalizing  the  re- 
sponses of  nonlinear  systems  for  a  class  of  inputs  to  the 
response  for  any  other  inputs.   The  response  to  one  type 
of  input  does  not  necessarily  contain  sufficient  informa- 
tion about  the  response  to  some  other  input.   Even  two 
inputs  which  are  of  the  same  waveshape  but  different  in 
amplitude  can  give  responses  that  are  entirely  different. 
This  makes  the  study  of  nonlinear  systems  quite  difficult. 
Even  with  this  difficulty,  there  is  little  choice  except 
to  tackle  the  problem  in  some  way  because  nonlinear  sys- 
tems occupy  a  very  important  role  in  practical  systems. 
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There  is  a  great  variety  of  possible  nonlinear  char- 
acteristics, and  there  are  different  ways  by  which  non- 
linear systems  are  classified.   One  classification  is  of  a 
mathematical  kind  based  on  the  equations  involved,  however 
the  study  of  nonlinear  systems  on  a  mathematical  basis  is 
extremely  difficult.   In  this  study,  types  of  nonlinearities 
encountered  in  servo  systems  are  considered,  namely  satura- 
tion and  hysteresis  effects.   The  consideration  of  nonlin- 
ear elements  in  the  operation  of  control  systems  have 
practical  significance  because  it  takes  linear  theory  nearer 
to  reality. 

B.    SUBHARMONIC  OSCILLATIONS 

A  linear  system  when  subjected  with  a  sinusoidal  in- 
put produces  a  steady-state  output  which  is  also  a  sinusoid 
and  oscillates  at  a  frequency  equal  to  the  input  frequency 
but  may  differ  from  the  input  in  amplitude  and  phase.   This 
is  true  for  all  frequencies.   However,  for  a  system  with 
suitable  nonlinear  characteristic  when  subjected  to  a  sin- 
usoidal input,  there  is  no  guarantee  that  the  steady-state 
output  is  also  a  sinusoid.   In  the  usual  case,  the  output 
wave  contains  frequency  components  which  are  integral  mul- 
tiples of  the  input  frequency.   Also,  the  system  can  pro- 
duce a  response  which  oscillates  at  a  frequency  which  is 
an  integral  submultiple  of  the  input  frequency.   This  fre- 
quency response  phenomenon  is  known  as  "subharmonic  reso- 
nance".  When  the  frequency  of  the  output  is  1/n  times  the 
driving  frequency,  where  n  is  an  integer,  then  the  output 
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oscillation  may  be  called  the  1/n  order  subharmonic  oscil- 
lation.  When  the  output  oscillates  with  the  input  fre- 
quency, the  output  oscillation  is  called  the  fundamental 
oscillation. 

The  subharmonic.  phenomenon  is  of  a  resonant  nature  and 
because  of  this  ,  the  subharmonic  component  is  often  of  large 
amplitude  which  causes  the  system  output  to  bear  little  re- 
semblance to  a  sinusoid  of  the  input  frequency.   The  output 
is  considerably  different  from  that  obtained  when  the  sub- 
harmonic phenomenon  does  not  occur.   For  systems  in  which 
the  output  is  intended  to  reproduce  the  input  as  closely  as 
possible,  as  in  the  motion  of  a  loudspeaker  cone  or  of  a 
servo  motor,  this  phenomenon  can  cause  a  distortion  that 
cannot  be  tolerated.   In  the  field  of  feedback  control  sys- 
tems, it  is  interesting  to  know  the  conditions  under  which 
subharmonic  oscillations  may  occur.   Subharmonic  oscilla- 
tions can  give  false  information  which  can  be  troublesome, 
so  it  is  desirable  to  eliminate  the  possibility  of  the 
occurrence  of  subharmonic  oscillations.   This  is  true  in 
the  case  of  recording  seromechanisms  which  are  frequently 
subjected  to  periodic  inputs.   The  possibility  of  the  occur- 
rence of  the  subharmonic  oscillations  reduces  the  reliability 
of  the  systems. 
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II.   DUAL- INPUT  DESCRIBING  FUNCTION 

The  requirement  that  a  nonlinearity  be  described  by 
an  equation  in  the  frequency  domain  and  that  this  equation 
be  compatible  for  use  with  the  transfer  function  of  the 
linear  components  gave  rise  to  approximating  the  effects 
of  the  nonlinear  component  by  using  a  linear  "approximate 
transfer  function"  or  describing  function.   The  approxima- 
tion is  made  by  defining  the  describing  function  in  terms 
of  the  Fourier  series  for  the  component  response  to  a  sin- 
usoidal input. 

The  describing  function  is  defined  as  the  ratio  of 
the  magnitude  of  the  fundamental  term  in  the  Fourier  se- 
ries for  the  output  wave  to  the  magnitude  of  the  input 
sinusoid  at  a  phase  angle  which  is  the  angle  between  the 
two  sine  waves  and  for  all  permissible  amplitude  and  fre- 
quency of  the  input  wave  [Ref.  9].   In  mathematical  form, 
let  the  input  be  ASinwt  and  the  fundamental  frequency  term 
in  the  Fourier  series  for  the  output  waveform  is 

F(A,w)Sin[cot+  (A,w)  ]  , 

then  the  describing  function  for  the  nonlinear  component 
is 


DF  =  lpCA^)  I    /<KA,u>) 
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Therefore,  the  describing  function  is  but  a  mathemati- 
cal linearization  of  a  nonlinearity  in  the  presence  of  a 
sinusoidal  input.   In  analysis  this  sinusoid  is  physically 
identified  as  the  system  limit  cycle  and  the  describing 
function  is  used  to  predict  its  amplitude  and  frequency. 
The  use  of  the  describing  function  is  based  on  the  assump- 
tion that  harmonics  generated  by  the  nonlinear  characteris- 
tics are  sufficiently  filtered  before  being  fed  as  input  to 
the  nonlinear  element. 

An  extension  of  the  describing  function  technique  is 
the  consideration  of  a  more  complex  signal  fed  to  the  non- 
linear element.   The  first  extension  is  the  use  of  dual- 
frequency  signal, 

v.  =  a  cos  (o)t+c|>)  +  b  cos  nut 

as  input  signal  to  the  nonlinearity  with  the  restriction 
that  a,  b,  and  n  must  be  real.   The  dual-input  describing 
function  can  be  defined  as  the  ratio  of  the  amplitude  of 
the  desired  frequency  component  in  the  output  waveform  to 
the  amplitude  of  the  component  of  the  same  frequency  in 
the  input  waveform.   Therefore,  in  the  determination  of 
the  dual-input  describing  function  it  is  necessary  to  find 
the  output  waveform  of  the  nonlinear  element  for  an  input 
composed  of  two  sinusoidal  signals  and  analyze  the  output 
to  find  the  components  of  a  particular  frequency.   In  the 
systems  considered  here,  these  signals  are  composed  of 
that  component  due  to  the  input  signal  at  a  fixed  frequency 
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and  the  component  due  to  the  fed  back  signal  which  may  be 
of  different  frequency  from  the  input  signal.   It  is  again 
assumed  that  the  harmonics  of  the  signal  fed  back  are  ade- 
quately filtered.   The  concept  of  the  describing  function 
is  thus  extended  to  allow  for  an  input  signal  composed  of 
two  sinusoidal  signals  of  different  amplitudes  and  fre- 
quencies . 

To  illustrate  the  derivation  of  the  dual-input  des- 
cribing function,  consider  the  cubic  characteristic, 


y  =  x3 


Let  the  input  be  of  the  form 

x   =   A  sin   ncot    +    B   sin    (cot  +  6) 

where 

A  and  B  are  independent  amplitudes , 

a)  and  nco  are  any  two  values  of  frequencies,  and 

6  is  an  independent  phase  angle. 

Substituting  the  value  of  x  in  the  nonlinear  equation,  the 
result  is, 

y  =  [A  sin  nwt+B  sin  (wt+0) ] 3 

=   A3sin3nojt  +  3A2Bsin2nwt   sin(wt  +  6) 

+  3AB2sin  ncot   sin2  (oot  +  9) +B3sin3  (wt  +  9) 
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Using  trigonometric  identities  to  express  the  above  ex- 
pression as  a  series  of  sine  waves  of  fundamental  and 
harmonic  frequencies  results  in 


y   =   A: 


■j  sin  ncot    -    j  sin    3nwt 


+  3A2B 


i  sin(cot  +  6)    -    j  sin{(2n+l)cot  +  0} 


+  j  sin{(2n-l)wt-6} 


+  3AB 


1  1 

y  sin   nwt    -    j  sin    {  (n+2)  oot  +  20 } 


j  sin    {(n-2)a)t-20) 


+   B 


3  1 

j  sin    (wt+0)    -    j  sin    3(cot  +  0) 


or 


•z  A  A  3 

y   =   —   (A2  +  2B2)    sin  noot    -    j-  sin   3nwt 


7T) 

+   ~    (B2  +  2A2)sin(oot  +  0)         ~-  sin    3(ojt  +  6) 


B_3 

4 


3A: 


sin{(2n-l)o)t-6}  -  sin{ (2n+l)wt+0 } 


3AB 


sin{  (n+2)wt+20>  +  sin{  (n-  2)  tot-  20} 
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By  inspection  of  the  above  expression,  the  sum  of  the  first 
and  third  terms  is  the  fundamental  input  function  equation, 
x  multiplied  by  the  factors  |-(A2  +  2B2)  and  |-(B2  +  2A2)  .   For 
the  other  terms,  there  are  values  of  n  which  have  to  be 
considered  as  will  be  explained  in  succeeding  paragraphs. 
Consider  the  following  cases: 

(1).   Non-harmonically  related  input  sinusoids.   In 
this  case  n  is  irrational,  i.e.,  n  is  not  expressible  as 
a  ratio  of  two  integers.   In  the  problem  being  considered, 
n  ^  1/3,  1,  or  3.   For  the  frequency  nu)/2iT,  the  dual- input 
describing  function  is  : 

|  A(A2  +  2B2)sin  neat 

N.(A,B)  =  _ . 

Av  '  J  A  sin  nut 

NA(A,B)  =  |  (A2+2B2) 

and  for  frequency  u)/2tt,  the  dual-input  describing  function 
is  : 

|  B(B2  +  2A2)sin(oot  +  9) 
VA,B)  =      B  sin(wt  +  e) 

NB(A,B)  =  |(B2+2A2) 

The  resulting  dual-input  describing  functions,  N. (A,B)  and 
Ng(A,B)  are  both  independent  of  the  phase  angle,  9  but 
dependent  on  the  input  signal  amplitudes  A  and  B. 

The  dual- input  describing  function  is  modified  if  n 
assumes  a  value  such  that  the  other  terms  in  the  series 
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have    the    frequency   of   either   of   the    input    signals    w/2tt   or 
nco/2ir,    i.e.,    n   =    1/3,    1,    or    3. 

(2).   Harmonically  related  input  sinusoids. 

(a)   n  =  3.   For  this  case  there  are  two  ad- 
ditional terms  in  the  output  waveform,  y  that  is  considered, 
i.e.,  the  fourth  and  eight  terms,  such  that  the  dual- input 
describing  function  becomes 

|^  (A2  +  2B2)  sin  3cot  -  ^-  sin(3wt  +  38) 
NaCA>B>3'6)  =  " A  sin  3u>t 

The  above  expression  shows  that  varied  phase  angles 
are  associated  with  the  terms,  so  the  phasor  representation 
is  utilized,  then 

f  A  (A2  +  2B2)eJ3Mt  -  |i  eH^t  +  3&) 
NA(A>B>3,6)  -  2 -33^ * 


Ae' 


NA(A,B,3,0)  =  |  (A2+2B2)  -  |^  ej30 


Similarly  for 


NB(A,B,3,9) 


S  (B2  +  2A2)e^wt  +  6^  -  3ABleJ(^-29) 


BeJCcot  +  0) 


NB(A,B,3,0)  =  |  (B2+2A2)  -  ^  e^39 

The  resulting  dual- input  describing  functions  are  complex 
and  dependent  on  A,  B,  and  6.   The  resulting  DIDFs  for  the 
other  two  cases,  n=l  and  n=l/3  as  shown  in  the  succeeding 
derivations  are  also  complex  and  dependent  on  A,  B,  and  6. 
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(b)   n=l.   For  this  case,  N.  and  NR  combine  to 
give  terms  of  frequency  oj/2tt.   Proceeding  in  the  same 
manner  as  for  the  case  n=3, 

|  A(A2  +  2B2)e^wt  +  |i  (B2  +  2A2)e^wt  +  9) 
NA(A,B,1,6)  = 


-A—-  Aeja)t 


3A2B   j(wt-e)    3AB2   j  (cot+26) 
,4     e 


Ae 


3  oat 


NA(A,B,1,6)  -  |  (A2+2B2)  +  ||  (B2+  2A2)e^6+  ^_  &-j 

+  3B^  eJ26 
4   e 


f  A(A2+2B2)ejwt  +  2  (B2  +  2A2)e^a3t  +  9) 

^,3,1,6)  "  " j(u,t+8) 

BeJ  v     } 


3A2B   j(ojt-e)    3AB2   j(ojt  +  26) 

4   e 4_  _^ 

,j(wt  +  e) 


Be' 


NB(A,B,1,6)  =  |4  (A2  +  2B2)  e"je  +  f(B2+2A2)  +  Ml  e^'2' 

+  MB   je 
4   e 


(c)   n  =  1/3.   This  case  is  identical  with  the 
case  for  n=3. 
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NA(A,B, 1/3,6)    = 


^    (A2+2B2)ej(wt/3)    -    3A2B    cJ  (tot/3  +  9) 


Ae 


j(o)t/3) 


NA(A,B,l/3,9)    =   |    (A2+2B2)    -    ^AB    eJ6 


NB(A,B,l/3,9)    = 


_AJ_  ejajt    +    3B    (B2  +  2A2)eJ(o)t  +  9) 


Be 


j  (cot+0) 


NB(A,B,l/3,    )    = 


A3      -19        3    r 
'    4B    e  +   T.    ( 


2p    (B2  +  2A2) 
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III.   SYSTEM  NONLINEARITIES 

In  most  practical  feedback  systems,  the  nonlinear 
element  is  in  the  forward  path.   The  basic  system  con- 
sidered in  this  study  is  shown  in  Figure  Ill-la,  where  N 
is  the  nonlinear  element  whose  gain  depends  on  the  magni- 
tude of  the  input  to  this  element  and  G(s)  is  the  linear 
frequency  dependent  portion. 

A.  SATURATION 

The  nonlinearities  considered  here  are  saturation, 
relay  with  hysteresis  and  backlash.   Most  physical  systems 
contain  certain  nonlinearities.   The  most  common  nonlinear 
phenomenon  is  characterized  by  a  saturating  or  limiting 
characteristic.   Limiting  is  almost  universally  present  in 
control  systems  since  most  instrumented  signals  can  take 
values  only  in  a  bounded  range.   Many  error  detectors, 
such  as  a  resolver  or  synchro  differential  have  a  restrict- 
ed range  of  linearity.   The  saturation  phenomenon  is  also 
present  in  amplifiers,  whether  electron  tubes,  transistors, 
magnetic  amplifiers  or  other  devices  are  used.   The  input- 
output  characteristic  is  shown  in  Figure  Ill-lb. 

B.  NON- IDEAL  RELAY 

The  other  nonlinearity  is  the  relay  with  hysteresis 
which  is  characterized  by  the  curve  shown  in  Figure  III-lc. 
The  hysteresis  delays  the  switching  operation  until  the 
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(a)   Block  diagram  of  basic  system. 
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(e)   Basic  system  with  two  nonlinear  elements 


FIGURE  III- 1 :   System  block  diagrams 
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output  is  past  the  desired  zero  point.   A  system  using  this 
type  of  relay  must  have  a  limit  cycle  in  steady  state. 

The  relay  is  a  device  to  apply  the  power  which  drives 
the  load.   In  a  relay  control  system,  the  power  amplifier 
is  a  relay  device.   The  relay  amplifier  is  desirable  be- 
cause it  can  be  simple,  rugged,  compact,  and  relatively 
cheap,  while  meeting  high  load  power  requirements.   Ex- 
amples where  relay  control  systems  are  typically  applied 
are  in  space  vehicle  attitude  control  and  aircraft  and 
missile  adaptive  control  systems. 

C.    BACKLASH 

Another  commonly  encountered  nonlinear  phenomenon  in 
most  mechanical  linkages  is  backlash.   In  control  systems 
these  are  present  in  most  mechanical  gears.   This  is  true 
even  for  new  ones  and  with  use,  the  backlash  increases. 
In  simple  physical  terms,  backlash  in  gears  is  the  free 
space  between  adjacent  teeth.   It  occurs  whenever  there 
is  inexact  linkage  between  mechanical  parts,  i.e.,  when 
there  is  slack  between  mating  gears.   The  backlash  char- 
acteristic is  shown  in  Figure  Ill-Id,  where  A  represents 
the  amount  of  backlash. 

In  one  way  or  another,  linear  and  nonlinear  systems 
will  encounter  the  backlash  problem  if  they  have  a  gear 
train.   In  the  field  of  control  systems,  a  system  is  in- 
cluded in  the  class  of  unstable  systems  if  it  continuously 
oscillates  with  constant  amplitude  and  zero  input  (limit 
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cycle).   Thaler  and  Pastel  [Ref.  9]  show  that  for  a  second 
order  servo  with  backlash,  the  existence  of  the  limit  cycle 
depends  on  the  value  of  the  damping  factor,  c .   For  z,    >   0.29, 
the  system  is  stable  so  that  for  £  <  0.29  the  effect  of 
backlash  is  to  cause  an  unstable  character.   If  the  phase 
portrait  is  drawn  for  the  system,  the  phase  portrait  will 
show  the  backlash  effect  near  the  origin.   For  systems 
where  it  is  desired  to  have  the  phase  trajectories  termin- 
ate at  the  origin,  the  backlash  effect  causes  serious 
trouble . 

The  three  nonlinearities  discussed  above  were  simula- 
ted on  the  electronic  analog  computer.   The  operation  of 
the  simulators  for  the  nonlinearities  are  explained  in 
Appendix  B.   The  characteristics  of  the  nonlinearities  were 
produced  very  nicely  as  shown  in  Figure  III- 2. 

Figures  III-3,  III-4,  and  III-5  show  the  block  dia- 
grams of  the  second  order  systems  with  the  indicated  non- 
linearity  which  were  simulated  on  the  analog  computer. 


24 


O) 


Backlash 


a 


m 


, 


m 


:.: 


:z.n 


i 


jxrr 


. 


tttt 


m 


m 


LLLUilULIIJ 


(c)   Relay  with  Hysteresis 


FIGURE  1 1 1 -  2  :   Simulation  output  for  nonlinearities 
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FIGURE  III-3a:   Block  diagram  of  second  order  system  with 

saturation. 


FIGURE  III-3b:   Analog  simulation  using  digital  logic  for 

simulation  of  saturation. 
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FIGURE  III-4a:   Second  order  system  with  non-ideal  relay 


FIGURE  III-4b:   Analog  simulation  of  system  with  relay 
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IV.   SUBHARMONIC  GENERATION 

Consider  the  basic  system  shown  in  Figure  Ill-la  with 
a  nonlinear  element  in  the  forward  path.   If  the  nonlinear 
element  is  replaced  by  an  element  of  unity  gain,  the  result- 
ing linear  system  is  assumed  to  be  stable. 

As  stated  earlier,  the  "Principle  of  Superposition" 
applies  to  linear  systems.   The  response  of  the  closed  loop 
system  to  a  given  signal  is  independent  of  any  other  signal 
which  may  be  applied  at  the  same  time.   When  a  nonlinear 
element  is  included  in  the  feedback  loop  this  is  no  longer 
the  case.   The  response  of  the  nonlinear  element  to  a  given 
signal  is  modified  if  another  signal  is  applied  at  the  same 
time.   The  effect  of  the  other  signal  may  be  considered  as 
varying  the  gain  of  the  nonlinear  element  for  the  primary 
input  signal.   The  case  of  interest  is  when  the  additional 
signal  input  is  also  sinusoidal  with  a  frequency  that  is 
an  integral  multiple  of  the  primary  input  frequency.   In 
In  this  case,  the  output  compcnent  of  the  angular  frequency, 
a)  may  be  modified  in  amplitude  and  phase  by  the  presence  of 
the  additional  input.   This  modifies  the  loop  gain  of  the 
system.   The  phase  change  is  caused  by  the  intermodulation 
products  generated  by  the  nonlinear  element. 

The  stability  of  the  closed  loop  system  is  affected  by 
varying  the  gain  of  the  nonlinear  element  since  the  nonlin- 
ear element  is  a  part  of  the  feedback  loop.   For  a  range  of 
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sinusoidal  inputs,  the  response  of  the  nonlinear  element 
may  be  so  modified  that  the  closed  loop  is  rendered  oscil- 
latory.  Where  the  system  is  such  that  the  open- loop  fre- 
quency response  plotted  on  the  complex  plane  do  not 
intersect  with  the  negative  real  axis,  the  nonlinear  ele- 
ments have  to  introduce  an  adequate  phase  change  in  order 
that  stable  oscillation  can  occur.   This  implies  that  the 
frequency  of  the  oscillation  must  be  related  to  the  input 
frequency.   What  is  of  interest  here  is  where  the  fre- 
quency of  the  oscillation  is  an  integral  submultiple  of 
the  primary  input  frequency  in  which  case  subharmonic  os- 
cillation is  generated. 
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V.   DISCUSSION 

The  systems  shown  in  Figures  1 1 1 -  3 ,  III-4,  and  III-5 
with  the  indicated  nonlinearities  were  simulated  to  in- 
vestigate the  occurrence  of  subharmonic  osciallations  and 
to  observe  the  effect  of  changing  the  variable  parameters 
in  the  systems  such  as,  the  damping  factor,  c  and  forward 

linear  gain,  K   in  Figures  1 1 1  -  3  and  1 1 1  -  5  and  the  amount 
&    >   m      & 

of  hysteresis,  2c  in  Figure  1 1 1  -  4 ,  in  the  generation  of 
subharmonic  oscillation. 

A.    SYSTEM  WITH  SATURATION  NONLINEARITY 

Figure  III-3a  shows  a  second  order  system  with  a  sat- 
uration nonlinearity .   The  nonlinear  characteristic  is 
defined  as : 

e=e.         -2<e.<2 

0      1  1 

e   =  2  e .  >  2 

o  1 

e   =  -2  e.  <  -2 

o  1 

and  the  gain  of  the  linear  portion  of  the  saturation  is 

unity.   The  transfer  function  of  the  linear  part  is, 

K 
G(s)  =  s(s  +  0.4) 

Considering  the  linear  part,  the  closed  loop  transfer 
function  is 

c  m 


Gn     sz  +  0.4s  +  K 
R  m 
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with  characteristic  equation 

s2  +  0.4s  +  K   =0 
m 

The  natural  frequency,  oo   and  the  damping  factor,  £  is 

determined  from  the  following  relationships: 


0)       =     vT~ 
n             m 

c  =  °*4  - 

0.2 

/K 
m 

In  the  simulation  of  the  system  in  the  electronic  analog 
computer,  K  was  varied  from  1  to  100  giving  a  variation 
in  £  from  0.2  to  0.02.   Typical  results  obtained  by  the 
simulation  are  shown  in  Figures  C-l  through  C-5. 

Ogata  [Ref.  6]  discussed  an  analysis  of  subharmonic 
responses  in  systems  with  a  nonlinear  element  whose  input- 
output  characteristic  curves  shows  saturating  or  limiting 
characteristics  and  showed  that  for  nonlinear  characteris- 
tic curves  that  are  skew  symmetric,  odd-order  subharmonic 
oscillations  can  occur. 

The  "28-degree  criterion"  for  predicting  responses  of 
saturating  feedback  systems  was  introduced  by  Douce  [Ref.  1] 
In  the  complex  plane,  two  lines  are  drawn  from  the  point 
s  =  - 1  on  each  side  of  the  negative  real  axis  and  each  line 
makes  an  angle  of  28°.   This  is  illustrated  in  Figure  V-l. 
The  region  bounded  by  the  two  lines  represent  the  critical 
region  for  a  system  with  saturation  nonlinearity .   Utilizing 
this,  the  Nyquist  diagram  of  the  linear  part  of  the  system 
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The  28-Criterion  (Douce  [Ref.  1]  .   Shaded  area 
shows  region  where  1/3  subharmonic  can  be  gen- 
erated in  systems  containing  a  limiter. 
Dotted  lines  represent  open  loop  frequency 
response  of  linear  system  with  indicated  value 
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provides  sufficient  information  to  predict  possible  occur- 
rence of  subharmonic  oscillations.   Subharmonic  resonance 
is  possible  when  the  open  loop  frequency  response  of  the 
linear  portion  enters  the  critical  region.   Figure  V-l 
shows  the  open  loop  frequency  response  of  the  linear  part 
for  different  values  of  K  .   Also  indicated  is  the  damping 
factor.   The  figure  shows  that  for  K   <  4.5  (s>  0.0944), 
the  open  loop  frequency  response  of  the  linear  portion 
lies  outside  the  critical  region  so  that  subharmonic  os- 
cillation cannot  occur.   For  K   >  4.5,  for  example  K   =50, 

m      *         r    m     ' 

the  linear  open  loop  frequency  response  is  inside  the  crit- 
ical region.   The  results  of  the  analog  computer  study 
agree  accurately  with  the  "28-degree  criterion"  as  can  be 
seen  in  Figure  V- 2  which  shows  the  domain  of  1/3  subhar- 
monic oscillation  as  a  function  of  input  frequency  and 
damping  factor,  £.   The  figure  also  shows  the  effect  of 
input  amplitude,  i.e.,  for  greater  input  amplitude,  the 
domain  decreases  and  the  system  has  to  be  more  lightly 
damped  to  sustain  subharmonic  oscillation. 

Subharmonic  oscillation  of  1/3  order  was  the  only  one 
observed  for  the  system  with  skew  symmetric  saturation  non- 
linearity.   Figure  V-3  shows  the  domain  of  1/3  subharmonic 
as  a  function  of  input  amplitude  and  input  frequency.   It 
can  be  seen  that  the  range  of  input  frequency  where  sub- 
harmonic oscillation  is  sustained  decreases  with  increasing 
input  amplitude.   The  same  observation  can  be  seen  in 
Figure  V- 2 .   It  can  be  deduced  that  subharmonic  oscillation 
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Input 

Amplitude 

"* 

■ 

6 

9 

16 

22 

30 

40 

Input 
Frequency 

4.4 

"■ 

• 

4.58 

80.0 

5.02 

75.0 

75.0 

5.65 

67.5 

68.75 

6.28 

57.5 

57.5 

55.0 

53.75 

6.91 

48.0 

46.5 

47.5 

47.5 

47.5 

7.53 

38.0 

40.0 

39.0 

40.0 

8.16 

34.0 

34.0 

32.5 

33.75 

8.79 

29.0 

29.0 

30.0 

30.0 

9.42 

25.0 

24.5 

27.5 

10.68 

20. 

0 

19.5 

18.75 

20.0 

10.99 

17. 

0 

17.0 

16.0 

17.5 

11.93 

16. 

0 

16.0 

15.5 

12.56 

14 

25 

13.75 

14.0 

13.19 

13 

0 

12.75 

12.5 

15.7 

9 

12 

9.0 

16.64 

7 

.75 

7.75 

18.84 

6 

.25 

20.1 

5 

.5 

21.35 

4 

.5 

21.66 

4 

.0 

Table  V-2  :   Output  amplitude  of  system  with  saturation 
C  =  0.02. 
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can  occur  only  for  a  range  of  input  amplitude.   For  a  small 
input  amplitude  such  that  the  error  signal  to  the  nonlinear 
element  is  less  than  the  saturation  level,  no  subharmonic 
oscillation  can  possibly  occur  because  the  system  will  just 
be  linear.   Also  shown  in  the  figure  is  the  variation  of 
output  amplitude.   For  a  particular  damping  factor,  the 
example  t,    =  0.02  in  Figure  V-4,  the  output  amplitude  varies 
as  shown  by  the  curve  irregardless  of  the  input  amplitude. 
The  output  amplitude  is  dependent  on  the  range  of  input 
frequency  that  the  subharmonic  oscillation  is  "locked  on". 
Figures  V-4  and  V-5  show  the  output  amplitude  variation 
for  different  values  of  £.   Ogata  [Ref.  6]  concluded  that 
subharmonic  oscillations  do  not  occur  as  response  of  the 
output  of  the  nonlinear  feedback  control  system  when  the 
output- input  amplitude  ratio  for  a  given  amplitude  of  a 
sinusoidal  input  is  2  ^  3  db  or  less.   This  is  equivalent 
to  a  ratio  of  1.256  ^  1.41.   The  results  of  the  simulation 
study  do  not  agree  with  this  criterion  as  can  be  seen  in 
Figures  V-4  and  V-5  where  the  lowest  output-input  amplitude 
ratios  are  1.25  and  1.125,  respectively.   The  possibility 
of  subharmonic  oscillation  is  greater  for  greater  output- 
input  ratio. 

Subharmonic  oscillation  of  1/2  order  was  observed  when 
the  saturation  nonlinearity  is  skew  asymmetric.   The  sat- 
uration characteristic  is 
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FIGURE  V-5 


Output  amplitude  for  1/3  subharmonic  due  to 
saturation.   Input  amplitude  =  16  volts. 
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The  domain  of  1/2  subharmonic  is  shown  in  Figure  V-6.   Fig- 
ures C-l  through  C-5  show  the  signal  waveforms  for  1/2  and 
1/3  orders  subharmonic  oscillations  and  also  the  waveforms 
of  the  error  signals.   Figure  C-2  shows  the  build  up  of 
1/2  order  subharmonic  and  Figure  C-4  shows  the  decay  of 
1/3  order  subharmonic  to  fundamental  oscillation. 

West,  Douce  and  Livesly  [Ref.  12]  presented  a  method 
of  using  the  dual- input,  describing  function  technique  to 
study  subharmonic  oscillation.   Figure  V-7  is  the  describ- 
ing function  loci  which  shows  the  response  of  a  saturation 
nonlinearity  to  a  sinusoidal  signal  of  amplitude  B  and  of 
a  signal  at  three  times  this  frequency  and  of  amplitude  A. 
Ref.  12  tabulates  other  describing  function  loci  for  dif- 
ferent values  of  amplitude  B. 

A  subharmonic  oscillation  of  1/n  order  where  n  is  an 
integer  will  oscillate  at  a  frequency  that  is  1/nth  of  the 
input  frequency.   To  investigate  this  subharmonic,  let  the 
input  to  the  nonlinearity  be 

A  cos  ncot  +  B  cos  (cot  +  9) 

The  signal  A  cos  nwt  is  assumed  to  be  known  and  the  des- 
cribing function  for  the  signal  B  cos  (oot  +  8)  is  superimposed 
on  the  open  loop  frequency  response  locus  of  the  linear 
part.   The  value  of  B,  such  that  the  describing  function 
cuts  the  frequency  response  locus  for  a  frequency,  w  is 
the  amplitude  that  corresponds  to  an  overall  open  loop  gain 
of  unity.   So  this  gives  a  possible  steady-state  amplitude 
of  the  resultant  subharmonic  oscillation. 
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The  amplitude  and  frequency  of  possible  1/3  order 
subharmonic  oscillation  is  found  as  follows.   Assume  that 
the  amplitude,  A  of  the  input  signal  is  known.   For  the 
system  to  oscillate  continuously  at  a  frequency  equal  to 
one- third  of  the  input  frequency,  the  overall  open  loop 
gain  must  be  unity.   Then  the  conditions  for  continuous 
oscillations  may  be  evaluated  by  means  of  the  describing 
function  technique.   The  input  to  the  nonlinear  element 
is 

A  cos  3wt  +  B  cos  (oat  +  e) 

Know  the  value  of  A,  the  describing  function  loci  that 
correspond  to  this  value  of  A  are  superimposed  on  the  open 
loop  frequency  response  locus  of  the  linear  system.   Refer 
ing  to  Figures  V-7a  and  V- 7b  this  shows  the  describing 
function  loci  for  A  =  4  and  A  =  8,  respectively.   The  open 
loop  frequency  response  locus  of  the  linear  portion  with 
a  transfer  function  of 

re    \    -       50 
G(;)W)  ~   juO+0.4) 

is  superimposed.   The  point  at  which  they  pass  through  the 
point  representing  the  frequency  response  at  the  subhar- 
monic frequency  determines  the  value  of  B,  which  is  the 
amplitude  of  the  subharmonic  oscillation. 

Figure  V- 8  shows  the  output  amplitude  for  1/3  sub- 
harmonic oscillation  for  a  second  order  system  with  sat- 
uration nonlinearity  where  the  saturation  level  is  ±  1. 
The  curves  were  obtained  from  the  results  of  the  analog 
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FIGURE  V-7a: 


Describing  function  loci  for  A  =  4. 
Shows  response  of  saturation  nonlin- 
earity  to  a  sinusoidal  signal  of 
amplitude  B  and  of  signal  AE  three 
times  this  frequency  and  of  ampli- 
tude A.  [Ref .  12] . 


46 


>Am 


50 


J^A 


oa -»-o.fi 


<^H\..^ 


FIGURE    V-7b 


Describing  function  loci  for 
A  =  8.   Shows  response  of  sat- 
uration nonlinearity  to  a  sin- 
usoidal signal  of  amplitude  B 
and  of  signal  at  three  times 
this  frequency  and  of  amplitude 
A.   Open-loop  frequency  response 
locus  of  linear  system  super- 
imposed. [Ref .  12 . ] 
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input 

Input  Fre-   A 
quency,  rad 

5.02 

5.65 

6.28 

6.91 

7.54 

8.16 

8.79 

9.42 
10.05 
10.68 
11.3 
11.93 
12.56 
13.19 
13.82 
14.13 
15.07 
15.7 
16.33 
16.96 
18.21 
18.84 


8  volts 


4  volts 


2  volts 


1  volt 


out 


err 


A 


out 


A 


err 


A 


out 


A 


err 


A 


out 


A 


err 


25.0   32.0 
19.0   25.0 

16.0   21.5   16.0  19.0 

12.5   18.0   12.5  15.5 

10.5   16.0   10.5  12.5   10.5   12.5 

9.0  11.0    9.0   10.5 


7.5 

9.5 

7.75 

9.0 

6.5 

8.5 

6.75 

8.0 

5.5 

8.0 

5.5 

7.5 

6.0 

6.5 

5.0 

7.5 

5.0 

6.0 

5.0 

5.5 

4.5 

5.5 

4.5 

5.0 

4.0 

5.0 

4.0 

4.5 

3.5 

4.5 

3.75 

4.1 

3.0 

4.25 

3.25 

4.0 

2.5 

4.0 

3.0 

3.5 

2.5 

3.75 

2.5 

2.35 

2.0 

2.0 

1.8 

1.6 

1.4 

3.0 

2.75 

2.5 

2.25 

2.25 

2.1 

2.0 

Table  V-3:   Amplitude  of  system  response  and  error  signal 
for  system  with  saturation  level  ±1. 
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computer  study.   Also  indicated  are  points  calculated  using 
the  describing  function  loci  of  Figure  V-7  and  those  tabu- 
lated in  Ref.  12. 

The  difficulty  with  the  dual- input  describing  function 
technique  is  the  derivation  of  the  dual-input  describing 
function  itself.   If  the  nonlinear  characteristic  can  be 
expressed  as  a  simple  power  series  or  simple  output-input 
relationship  such  as  the  cubic  nonlinearity  considered  in 
Chapter  II,  the  dual-input  describing  function  can  be  de- 
termined by  simple  algebraic  analysis.   However,  in  the 
case  of  practical  feedback  systems  where  the  nonlinearity 
is  rarely  expressible  as  a  simple  power  series,  like  in  the 
case  of  saturation,  relay,  and  backlash,  the  algebraic  an- 
alysis of  the  nonlinear  element  becomes  difficult.   Graph- 
ical methods  are  possible  but  very  laborious.   Due  to  the 
complexity  of  these  nonlinearities ,  the  use  of  computers 
for  determining  the  relevant  dual-input  describing  func- 
tion is  necessary. 

For  the  purpose  of  this  study,  the  derivation  of  the 
dual-input  describing  function  for  the  nonlinearities  are 
not  included.   The  DIDF  derivation  can  be  treated  as  a 
separate  problem  and  the  results  can  then  be  applied  to 
any  form  of  linear  systems  to  study  the  effect  of  the 
nonlinearity. 

B.    SYSTEM  WITH  NON- IDEAL  RELAY 

The  other  second  order  system  simulated  to  observe 
the  occurrence  of  subharmonic  oscillation  contained  a 
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double  valued  nonlinearity .   Figure  III-4a  shows  a  second 
order  system  with  a  relay  with  hysteresis.   Referring  to 
this  figure,  the  nonlinearity  can  be  defined  by  the  fol- 
lowing equations.   Assume  the  input  to  the  nonlinearity 
is 


e.  =  A  sinojt 


then 


e 

0 

-a 

0 

<   tot   <   3 

e 

0 

a 

3 

<     U)t     <     7T  +  3 

where 


3  =  sin   c 


A 
The  conventional  describing  function  [Ref.  9]  is 


4a   /   .  -1 
^A  /  Sin  £ 


and  the  plot  of  the  inverse  of  this  describing  function 
for  a  =  1  (constant)  and  three  different  values  of  c,  where 
2c  is  the  amount  of  hysteresis,  is  shown  in  Figure  V-9. 
Superimposed  on  the  -1/G^  curve  is  the  plot  for  G(jco)  of 
the  linear  portion  of  the  system.   The  figure  shows  that 
for  a  particular  amount  of  hysteresis,  the  point  of  inter- 
section between  the  G(jw)  curve  and  the  relevant  -1/G~ 
curve  represents  a  limit  cycle  whose  amplitude  and  fre- 
quency can  be  determined  easily.   For  example,  consider 
c  =  0.2,  the  limit  cycle  frequency  is  approximately  1.8 
rad/sec  with  an  amplitude  of  0.46.   The  Bode  plot  of  the 
simple  linear  portion  is  shown  in  Figure  V-10. 
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FIGURE  V-9:   Describing  function  for  relay  with  hysteresis 


52 


o 

CO 

i 


O 
I 


=8  $ 


o 


v / 

u 

m 
o 

o 

rH 
Ph 

(1) 

o 


i 
> 

W 
Pi 

1—1 

PL, 


o 
on! 


o 

ID 


o 

I 


o 

I 


53 


The  system  is  Figure  III-4a  has  a  sinusoidal  input  of 
known  amplitude  and  frequency.   The  output  of  the  relay, 
f (t)  is  assumed  to  lose  its  higher  frequency  components  due 
to  the  low  pass  characteristics  of  G(s) .   Then,  the  input 
to  the  nonlinearity  can  be  written  as 

A  cos  nwt  +  B  cos  (wt+9) 

where  A  cos  nwt  is  the  primary  input  signal  of  known  am- 
plitude and  frequency.   The  dual- input  describing  function 
of  subharmonic  oscillation  developed  by  West,  Douce  and 
Livesly  [Ref.  12]  can  be  used  if  the  dual-input  describing 
function  is  available. 

In  the  analog  computer  study  for  the  system,  subhar- 
monic oscillation  as  high  as  1/2  order  and  as  low  as  1/13 
order  were  observed.   The  subharmonic  oscillation  once  ini- 
tiated is  stable  and  stays  "locked  on"  a  reasonable  range 
of  frequency.   There  is  an  overlapping  of  domains  of  dif- 
ferent subharmonic  order,  as  shown  in  Figure  V-ll.   The 
subharmonic  oscillation  that  is  stable  is  the  subharmonic 
order  initiated  first.   It  can  be  seen  from  Figure  V-ll 
that  odd-order  subharmonics  have  a  much  larger  domain  than 
even-order  subharmonics  and  that  the  domain  decreases  for 
lower  odd-order  subharmonics.   The  figure  also  shows  the 
amplitude  of  the  output  which  is  the  amplitude  of  the  sub- 
harmonic oscillation.   This  is  the  amplitude  curve  for  all 
subharmonic  modes.   If  the  output  oscillates  at  a  frequency 
of  2.45  rad/sec,  then  the  output  amplitude  is  0.225  irre- 
gardless  of  the  subharmonic  mode. 
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Input 

Ampli 

i 

0.1 

0 

.12 

0 

.18 

1 

3.21 

0 

.27 

tude 

Order 

to- 
m 

A 
out 

in 

A 
out 

0)  . 

in 

A 
out 

O). 

in 

A 
out 

O). 

in 

A 
out 

1/2 

3.77 

.45 

3.77 

.40 

3.77 

.375 

3.77 

.375 

■3.77 

.375 

4.08 

.40 

4.08 

.35 

4.08 

.338 

1/3 

5.02 

.538 

5.02 

.543 

5.02 

.52 

5.02 

.55 

5.02 

.513 

5.65 

.40 

5.65 

.40 

5.65 

.4 

5.65 

.40 

5.65 

.40 

6.28 

.338 

6.28 

.35 

6.28 

.325 

6.28 

.325 

6.28 

.325 

6.91 

.25 

6.91 

.275 

6.91 

.26 

7.22 

.225 

7.54 
8.15 

.213 
.175 

7.54 

8.16 

8.79 

9.42 

10.05 

10.68 

11.3 

11.93 

.213 

.175 

.15 

.125 

.113 

.10 

.088 

.075 

1/4 


1/5 


1/6 


1/7 


6 

.91 

.475 

6, 

.91 

.475 

7, 

.54 

.4 

8 

.05 

.363 

8.16 

.325 

7, 

.53 

.40 

7 

.53 

.4 

7 

.85 

.35 

8, 

.16 

.538 

8, 

.16 

.525 

8 

.16 

.538 

8 

.16 

.55 

8.79 

.45 

8, 

.79 

.45 

8, 

.79 

.45 

8, 

.79 

.45 

9, 

,42 

.375 

9.42 

.375 

9. 

,42 

.375 

9, 

,42 

.375 

9, 

.42 

.388 

10, 

,05 

.31 

10.05 

.33 

10, 

,05 

.325 

10, 

.05 

.338 

10, 

.68 

.28 
.26 
.23 
.2 

11.3 

11.93 

12.56 

13.19 

13.82 

14.44 

15.07 

15.7 

.26  5 

.24 

.2 

.188 

.175 

.163 

.15 

.125 

10. 

,05 

.5 

10, 

,68 

.425 

11, 

.3 

.4 

11, 

,62 

.375 

10. 

,67 

.44 

11, 

,3 

.375 

11, 

.93 

.34 

11, 

.93 

.34 

11, 

,3 

.38 

12, 

,56 

.31 

11. 

,93 

.49 

11, 

,93 

.49 

11, 

,93 

.49 

12, 

,87 

.4 

15.7 

.26 

12. 

,56 

.44 

12, 

,56 

.44 

12, 

,56 

.43 

13, 

,82 

.35 

16.33 

.25 

13. 

,19 

.39 

13. 

,19 

.39 

13, 

.19 

.39 

14, 

,44 

.33 

16.96 

.225 

13. 

,82 

.35 

13, 

,82 

.35 

13, 

,82 

.35 

15, 

,07 

.29 

17.58 

.21 

14, 

.44 

.31 

15, 

,7 

.26 

18.21 

.2 

15, 

.07 

.3 

16, 

,33 

.25 

18.84 

.19 

16, 

,01 

.25 

16, 

,96 

.24 

20.1 
21.35 

.15 
.14 

Table  V-4:   Simulation  Data  for  Non-Ideal  Relay; 
a  =  1 ,  c  =  0 . 2 
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INPUT        n 
AMPLITUDE 

1 

0.12 

0.18 

0.21 

0.27 

Order 

a).   -A 
in 

out 

to.   A   . 
in   out 

a).   A 
in   out 

03.    A 

m   out 

a).   A 
in   out 

1/8 

14.44 

.43 

14.44  .43 

15.7   .35 

16.01 

.35 

15.07 

.39 

15.0  7  .39 
15.7   .35 

16.33  .33 
16.96  .3 
17.58  .28 

17.27 
18.21 
18.84 

.28 
.26 

.24 

1/9 

15.7 

.46 

16.33  .43 

16.96  .39 

18.21 

.34 

16.33 

.43 

16.96  .4 

17.58  .36 

19.15 

.3 

16.96 

.38 

17.58  .36 

18.84  .31 

20.1 

.26 

17.58 

.35 

18.21  .35 

19.47  .29 
20.41  .25 

20.72 
21.35 

.25 
.23 

1/10 

18.21 
18.84 
19.47 

.43 
.39 
.36 

18.84  .39 

1/11 

20.09 
20.72 
21.35 

.44 
.39 
.36 

INPUT 

AMPLITUDE    0. 

3 

0.375 

0.4 

0.45 

0.51 

Order 

1/3 

4.7 

.72 

5.02  .58 

4.71  .63 

4.71 

.63 

4 

.9   .53 

5.02 

.55 

5.65  .4 

5.02  .53 

5.02 

.55 

5 

.15  .45 

5.65 

.43 

6.28  .33 

5.65  .4 

5.65 

.38 

5 

.27  .44 

6.28 

.33 

6.91  .28 

6.28  .33 

6.28 

.33 

5 

.4   .41 

7.54 

.21 

7.54  .21 

6.91  .25 

5 

.53  .39 

8.79 

.15 

8.16  .17 

7.54  .2 

10.05 

.11 

11.3 

.09 

. 

12.56 

.07 

13.19 

.06 

1/5 

9.42 
10.05 
10.68 
13.82 
15.07 
16.33 
16.96 

.38 

.34 

.3 

.18 

.14 

.12 

.11 

Table  V-5:   Simulation  Data  for  Non-Ideal  Relay  (continuation) 
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.1  ■■ 


G(s)  = 


S*  +  S  +  1 


Boundary  points  obtained 
by  analog  computer. 


z0— 


V 

A 


1/3  Subharmonic 

1/5  Subharmonic 


/•  -i  \ 

I    /  *•'/  y    1/7  Subharmonic 

/  •  v 


■\  Output  Amplitude 


2         3         4 
OUTPUT  FREQUENCY,  rad/sec 


5" 


FIGURE  V-ll :   Domains  for  subharmonic  output  for  system 
with  relay:  a  =  1,  c  =  0.2. 
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Figure  V-12  shows  the  domain  of  subharmonic  oscilla- 
tions for  a  =  1  and  c  =  0.1  adapted  from  the  work  of  Sawaragi 
and  Akashi  [Ref .  8] .   The  points  plotted  represent  data  that 
were  obtained  from  the  analog  computer  experiment. 

Figures  C-6  through  C-16  show  the  simulation  results 
for  the  system  with  non-ideal  relay.   In  the  system  with 
saturation  nonlinearity ,  subharmonic  generation  was  not 
self  starting.   The  occurrence  of  subharmonic  oscillation 
was  dependent  on  the  initial  conditions,  such  as,  the  sud- 
den change  of  input  frequency  or  input  amplitude.   For  the 
system  with  non-ideal  relay  nonlinearity,  subharmonic  oscil- 
lation was  self  starting.   This  can  be  seen  in  Figure  C-6 
which  shows  the  system  response  changing  from  fundamental 
oscillation  to  1/3  subharmonic  oscillation.   Figure  C-6 
also  shows  the  relationship  between  the  amplitude  of  the 
output  when  oscillating  at  the  fundamental  frequency  and 
when  oscillating  at  the  subharmonic  frequency.   Figure  C-10 
shows  the  change  from  subharmonic  to  fundamental  oscilla- 
tion, i.e.,  from  1/5  order  subharmonic  to  fundamental.   The 
amplitude  of  the  output  at  the  fundamental  frequency  is 
practically  negligible  compared  with  the  output  amplitude 
when  oscillating  at  a  subharmonic  frequency.   Figures  C-7 
through  C-16  show  the  changing  modes  of  the  subharmonic  os-' 
dilation  due  to  change  in  input  frequency.   There  were  in- 
stances when  the  output  oscillates  at  two  different  subharmonic 
modes  at  the  same  input  frequency.   Figure  C-16  shows  the 
output  oscillating  at  two  different  subharmonic  modes;  name- 
ly, 1/11  order  and  1/13  order  at  the  same  input  frequency. 
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C.    SYSTEM  WITH  TWO  NONLINEARITIES 

The  third  system  simulated  to  investigate  the  occur- 
rence of  subharmonic  oscillation  contained  saturation  and 
friction  controlled  backlash  nonlinearities  as  shown  in 
Figure  III-5a.   It  is  a  two  loop  system  with  a  tachometer 
feedback  for  the  inner  loop  to  provide  the  necessary  damp- 
ing to  make  the  system  stable. 

Using  standard  procedure,  the  characteristic  equation 

is 

s2  +  K  K^s  +  K   =0 
m  t     m 

It  is  clear  that  the  damping  factor,  £  is  dependent  on  K 


in 


and  K  ,  i.e., 


Kt 


s    2    m 

Han  [Ref.  4]  presented  a  procedure  to  predict  the 
existence  of  limit  cycles  in  a  system  with  two  nonlinearities 
such  as  the  system  shown  in  Figure  II  I- 5a.   For  the  purpose 
of  observing  the  occurrence  of  subharmonic  oscillation  in 
the  system,  assume  that  the  system  oscillates  with  a  limit 
cycle  so  that  the  input  to  the  first  nonlinearity  is  of  the 
form 

A  cos  nwt  +  B  cos  (wt+0) 
where  A  cos  nwt  is  the  primary  input  signal  of  known  ampli- 
tude and  frequency.   The  linear  frequency  dependent  ele- 
ment, G(s)  is  a  low-pass  filter  which  is  assumed  to  adequately 
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attenuate  the  higher  frequencies  so  that  harmonics  can  be 
neglected  in  e.  .   It  also  follows  that  the  input  to  the 
second  nonlinearity  is  a  sinusoid. 

The  analog  computer  simulation  was  performed  for  dif- 
ferent values  of  Km  and  K   and  the  results  are  tabulated 
in  Table  V-6  and  V-7.   Figures  C-17  through  C-23  show  the 
signal  waveforms  obtained.   Figure  C-17  and  C-18  show  the 
waveforms  for  1/3  and  1/5  orders,  respectively.   It  can 
be  seen  that  the  signal  e.   which  is  the  input  to  the  back- 
lash nonlinearity  oscillates  at  the  subharmonic  frequency, 
i.e.,   in  Figure  C-17,  e.   is  oscillating  at  a  frequency 
that  is  1/3  of  the  input  frequency  and  in  Figure  C-18,  at 
1/5  of  the  input  frequency.   The  system  output,  6   oscil- 
lates at  the  same  frequency  as  e.   with  a  phase  difference 
of  180°.   This  shows  that  the  subharmonic  oscillation  is 
due  to  the  effect  of  the  saturation  nonlinearity.   The  ef- 
fect of  the  backlash  nonlinearity  can  be  analyzed  by  using 
the  conventional  describing  function  since  the  input  to  the 
backlash  nonlinearity,  e.   is  only  composed  of  one  sinusoid- 
al signal.   Figure  C-19  shows  the  change  of  fundamental  os- 
cillation to  1/3  subharmonic  oscillation  and  Figures  C-20 
and  C-21  show  the  build  up  of  subharmonic  oscillations. 
It  can  be  seen  also  from  these  figures  that  the  subhar- 
monic oscillations  are  self  starting  as  compared  with  the 
system  with  saturation  nonlinearity  only. 

The  domains  of  subharmonic  oscillation  as  a  result  of 
the  simulation  study  are  shown  in  Figures  V-13,  V-14,  and 
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V-15.   The  domain  of  the  1/3  subharmonic  oscillation  is  of 
the  order  of  the  natural  frequency.   For  lower  order  sub- 
harmonic  modes,  i.e.,  1/4,  1/5,  1/6,  1/7,  ...,  subharmonic 
oscillation  occurs  at  a  higher  input  frequency. 
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Input 

Amplitude     3.0         9.0       21.0        30.0 

M     Out-       Out-       Out-         Out- 
Order     in   put,v    in  put ,v    in  put,v    in   put,v 

1/2  20.41  4.75  11.93  15.0 
20.72  4.75  12.56  14.0 
21.35  4.5  13.82  13.0 
25.12       15.07  12.0 

16.33  12.0 

17.58  11.75 

1/3   22.61  15.5  11.3  66.7  7.54  135.0   8.16   130. 

23.86  13.0  12.56  60.  8.16  125.0   8.79   125. 

25.12  11.75  13.82  48.  8.79  105.    9.42    90. 

26.38  10.  15.02  39.  9.42   85.   10.68    72.5 

27.63  8.88  16.33  32.  10.05   80.   11.93    60.0 

28.89  8.25  17.58  27.  10.68   75.   13.19    50. 

30.14  7.5  18.84  22.  11.93   57.5  13.82    42.0 

31.4  6.1  20.1  20.  12.56   50.   15.07    37.5 

36.66  6.0  21.35  17.5  13.82   45. 

33.91  5.5  22.61  15.0  15.07   37.5 

35.17  4.5  23.86  13.0  16.33   30.0 

28.26  9.5  17.58   25.0 

29.52  8.0  18.84   20.0 

1/4   35.8         30.14  10.0  18.12        15.7  27.5 

31.4    9.0  18.79       16.33  27. 

34.54   8.5  20.1         16.96  25.0 

21.98       17.58  23. 

18.21  22 

1/5  37.68   9.0  22.61        11.93 

.  13.19 

t0         14.44  117.5 

28.26       15.7  97.5 

16.96  80 

18.21  72.5 

20.1  62.5 

21.35  52.5 

22.61  47.5 


Table  V-6:   Domain  of  subharmonic  for  system  with  two 
nonlinearities ,  c  =  0.283,  oon  =  14.14 
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VI.   CONCLUSIONS 



The  investigation  of  subharmonic  oscillation  in  non- 
linear systems  was  undertaken  by  means  of  an  analog  com- 
puter study  on  three  nonlinear  second  order  feedback 
control  systems.   The  following  conclusions  are  summarized 
from  the  result  of  this  study. 

When  the  output  of  the  system  with  saturation  non- 
linearity  is  oscillating  with  the  input  frequency,  sub- 
harmonic  oscillation  cannot  occur  unless  some  kind  of  sudden 
disturbance  is  given  to  the  input.   This  sudden  disturbance 
may  be  a  sudden  disturbance  may  be  a  sudden  increase  in 
input  amplitude  or  frequency.   This  is  not  the  case  for 
the  systems  with  hysteresis  nonlinearity .   In  the  system 
with  non-ideal  relay  and  backlash  characteristics,  sub- 
harmonic  generation  was  self  starting. 

The  possibility  of  the  occurrence  of  subharmonic  os- 
cillation is  great  when  the  output-input  amplitude  ratio 
is  large.   In  the  case  of  saturation  nonlinearity,  there 
is  a  minimum  output- input  amplitude  ratio  for  subharmonic 
oscillation  to  occur,  i.e.,  the  output-input  amplitude 
ratio  is  greater  or  equal  to  1.125.   For  the  system  with 
non-ideal  relay,  subharmonic  oscillation  can  occur  for 
output-input  amplitude  ratio  less  than  unity. 

The  "28-degree  criterion"  introduced  by  Douce  [Ref.  1] 
gives  accurate  results  in  predicting  the  possibility  of 
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subharmonic  responses  in  a  system  with  saturation  nonlin- 
earity.   This  is  a  useful  technique  in  the  elimination  of 
subharmonic  oscillation  for  systems  with  a  limiting  non- 
linearity  . 

Subharmonic  oscillation  once  initiated  is  sustained 
over  a  reasonable  range  of  frequency.   There  is  an  over- 
lapping of  these  frequency  ranges  for  different  subhar- 
monic modes.   In  the  region  where  two  or  more  orders  of 
subharmonic  oscillation  can  occur,  the  order  that  occurs 
first  dominates  and  tends  to  continue  and  is  quite  stable. 
The  frequency  of  the  subharmonic  oscillation  has  a  minimum 
depending  on  its  order.   The  frequency  for  1/3  order  sub- 
harmonic is  usually  in  the  neighborhood  of  the  natural 
frequency  of  the  system.   Subharmonic  oscillation  is  de- 
pendent on  the  input  signal  amplitude.   It  can  occur  only 
within  a  certain  range  of  signal  amplitude.   It  can  be  con- 
cluded that  subharmonic  oscillation  can  occur  only  with  a 
certain  range  of  input  amplitude  and  frequency. 

The  system  with  a  limiting  nonlinearity  is  dependent 
on  the  damping  factor  for  subharmonic  oscillation  to  occur. 
The  system  has  to  be  very  lightly  damped,  i.e.,  z,  <  0.1  in 
order  for  the  system  response  to  oscillate  at  a  subharmonic 
frequency.  If  the  input  amplitude  is  increased,  the  system 
has  to  be  much  more  lightly  damped  for  subharmonic  oscilla- 
tion to  occur. 

The  dual- input  describing  function  technique  developed 
by  West,  Douce,  and  Livesly  [Ref.  12]  appears  to  be  a  very 
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useful  technique  for  predicting  quantitatively  the  ampli- 
tude and  frequency  of  subharmonic  oscillation.   Recommend- 
ed for  future  study  is  the  derivation  of  relevant  dual- input 
describing  functions  (two  sinusoidal  input)  and  application 
of  the  results  to  any  linear  system  with  the  relevant  non- 
linearity.   Quantitative  analysis  of  the  occurrence  of  sub- 
harmonic  oscillations  should  be  undertaken. 
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APPENDIX  A:   PROGRAM  SYMBOLS 


(1)   Servo  Set  Potentiometer 


e . 

1 


K 


e   =  Ke .  where  K  <  1 
o      1 


(2)   Summer  Amplifier 


e   =  -  (e..  +  10e.o) 

o      *■  ll      \2J 


(3)   Inverter 


(4)   Integrator 


e   =  -  e . 
o       1 


e   =  -    (e. , +e. ~  +  10e . _)dt 
o       v  ll  i2    ±y 


(5)   Comparator 


X 

Y 


C 
C 


Output  Logic: 

Algebraic  sum  of  X  +  Y  >  0 

C  =  -6  volts  (true) 

C  =   0  volt   (false) 
Algebraic  sum  of  X  +  Y  <  0 
C  =   0  volt   (false) 
C  =  -6  volts  (true) 
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(6)   Digital/Analog  (D/A)  Switch 


Logic  Control: 

The  switch  is  closed  if  either 
-6  volts  is  applied  or  the 
input  is  left  open;  the  switch 
is  open  if  0  volt  is  applied. 


(7)   OR  Gate 


+ 


-4  +  L2 


(8)   AND  Gate 


■Ll  +  L2 
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APPENDIX  B:   ANALOG  SIMULATION 

1.    Nonlinearities 
(a)   Saturation 

Figure  B-la  shows  the  saturation  characteristic 
curve  and  Figure  B-lb  is  the  circuit  diagram  for  the  simu- 
lation using  the  analog  computer.   Referring  to  Fig.  B-lb, 
the  simulation  works  as  follows : 

(1)  e.  <  -  a 

v  J         1 

Algebraic  sum  of  CI  <  0  -*■  CI  false,  CI 
true  and  also  the  algebraic  sum  of  C2  <  0  ■*■  C2  false, 
C2"  true  such  that  the  output  of  the  AND  gate  =  0  and 
D/Al  is  open.   For  CT  true,  D/A2  is  closed  and  for  C2 
false,  D/A3  is  open,  therefore, 

e   =  +  V, 
o      1 

(2)  Algebraic  sum  of  CI  >  0  ->  CI  true,  CT 
false  and  the  algebraic  sum  of  C2  <  0  -*■   C2  false  and  C2 
true.   The  output  of  the  AND  Gate  =  -6  volts  and  D/Al 
is  closed.   For  CT  false,  D/A2  is  open  and  C2  false, 
D/A3  is  open,  therefore, 


e   =  -  me . 
o        1 


(3)   0  <  ei  <  b 

Algebraic  sum  of  CI  >  0  -*■  CI  true,  CI 
false  and  the  algebraic  sum  of  C2  <  0  -*■  C2  false  and  C2 
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Vt   , slopes 


VYL 


-*c? 


FIGURE  B-la:   Saturation  characteristic  curve. 


*-  O/Al 

*»  D/A5 


b/ioo 


FIGURE   B-lb: 


Analog  simulation  for  saturation.   The 
symbol  "g"  refers  to  frid  of  amplifier 
and  100  is  the  reference  voltage  of  the 
computer  being  used. 
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true.   The  output  of  the  AND  Gate  =  -6  volts  and  D/Al  is 
closed.   For  CT  false,  D/A2  is  open  and  C2  false,  D/A3 
is  open,  therefore, 

e   =  -  me . 
o       1 

(4)   e.  >  b 

v  J        i 

Algebraic  sum  of  CI  >  0  ■*  CI  true,  CT 
false  and  the  algebraic  sum  of  C2  >  0  -*■  C2  true  and  C2" 
false.   The  output  of  the  AND  Gate  =  0  and  D/Al  is  open. 
For  CT  false,  D/A2  is  open  and  for  C2  true,  D/A3  is 
closed,  therefore, 

e   =  -  V0 
o      2 

In  the  description  of  operation  above,  -6  volts  represents 
a  logical  1  and  0  volts  represents  a  logical  0.  The  logic 
control  for  the  D/A  switches  is  the  output  of  the  compara- 
tor or  of  the  AND  Gate. 

(b)   Relay  with  Hysteresis 

Figure  B-2a  is  the  characteristic  curve  of  a 
non- ideal  relay  with  hysteresis  and  Figure  B-2b  is  the 
diagram  for  analog  simulation.   The  Track-Track  (T/T)  unit 
is  composed  of  two  D/A  switches  with  a  common  G  terminal. 
Figure  B-2c  is  the  analog  circuit  for  the  T/T  unit.   Switch 
T  is  closed  and  T  is  open  when  a  logical  1  is  applied  to 
the  input  or  the  input  is  left  open.   When  the  input  is  a 
logical  0,  the  switch  states  are  reversed.   When  the 
switches  are  open,  the  resistor  switch  ends  are  grounded. 
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FIGURE  B-2a:   Characteristic  curve  of  relay  with  hysteresis 


FIGURE  B-2b:   Analog  simulation  of  relay  with  hysteresis 
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AAAA S~ 
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FIGURE  B-2c:   Analog  circuit  of  track-track  (T/T)  unit 
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(c)   Backlash 

The  backlash  characteristic  curve  is  shown  in 
Figure  B-3a  and  the  analog  simulation  circuit  diagram  is 
shown  in  Figure  B-3b.   The  backlash  characteristic  is 
achieved  as  follows: 

As  e.  increases  from  zero  nothing  happens  to  e 
until  (e .  -  e  )  =  e'  >  a.   At  that  point  the  integrator 
output  will  begin  to  change  at  the  rate  m(e.  -  e  ) ,  and 
e   will  approach  (e.  -  a)  through  the  feedback  correction 
of  the  loop.   As  soon  as  e.  changes  direction,  (e-  -  e  ) 

r  1  b  '   ^  1      O 

becomes  less  than  a  resulting  in  e"  going  to  zero  and  the 
output  of  the  integrator  holding  its  value.   As  e.  con- 
tinues to  decrease,  e   remains  unchanged  until  e'  becomes 

'   o 

less  than  b,  at  which  time,  the  output  of  the  integrator 
begins  to  change  at  the  rate  m(e.  -  e  ).   Again  when  e. 
changes  direction,  e'1  goes  to  zero  and  the  output  of  the 
integrator  remains  at  some  steady  voltage. 

The  results  of  the  simulation  for  the  three  nonlin- 
earities  considered  are  shown  in  Figure  III-2. 
2.   Scaling 

When  simulating  problems  on  an  analog  computer,  scalin] 
is  quite  important.   There  are  many  approaches  to  this 
technique.   What  is  of  interest  here  is  to  point  out  that 
it  is  necessary  to  perform  magnitude  and  frequency  scal- 
ing on  the  analog  simulation  of  the  nonlinearities  simula- 
ted in  this  study.   This  is  true  for  the  two  nonlinearities 
with  hysteresis.   The  simulation  for  the  saturation  nonlin- 
earity  only  requires  amplitude  scaling. 
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FIGURE  B-3a:   Backlash  characteristic  curve 


FIGURE  B-3b:   Analog  simulation  for  backlash 
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FIGURE  C-l:   1/2  subharmonic  oscillation  due  to  asymmetric 

saturation  input  =  6  Sin  1.6irt:  K   =  20 

r  m 
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FIGURE  C-2:   Fundamental  to  1/2  subharmonic  due  to  asymmetric 
saturation,  smm/sec. 
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FIGURE  C-3:   1/3  subharmonic  oscillation  due  to  symmetrical 

saturation  input  =  3  Sin  5.67Tt,  K   =  100 

r  m 
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FIGURE  C-4: 


1/3  subharmonic  oscillation  to  fundamental  os 
dilation  input  frequency  from  2.9hz  to  3.0hz 


K   =  70 
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FIGURE  C-5:   1/3  subharmonic  oscillation  due  to  saturati 

input  =  3  Sin  5.2irt,  K   =  70 
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FIGURE  C-8:  1/3  to  1/5  order  subharmonic  for  system  with 
relay;  a  =  1,  c  =  0.1.  Input  frequency  from 
1.8  hz  to  1.9  hz. 
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FIGURE  C-10:   1/5  subharmonic  oscillation  to  fundamental 
oscillation.   a  =  1,  c  -.01.   Input  fre- 
quency from  4.0  hz  to  4.6hz. 
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FIGURE  C-ll:   1/5  to  1/7  subharmonic  oscillation  for  relay 
a  =  1,  c  =  0.1.   Input  frequency  from  3.1hz 
to  3 . 2hz . 
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FIGURE  C-12:   1/7  to  1/5  subharmonic  oscillation.   Input 

frequency  from  2.5hz  to  2.4hz,  a  =  1,  c  =  0.1 
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FIGURE  C-13:   Change  of  input  frequency  from  3.2hz  to  3.4hz 
causing  change  of  subharmonic  oscillation  from 
1/7  to  1/9;  a  =  1,  c  =  0.1. 
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FIGURE    C-14 


1/9    to    1/7    subharmonic   for   relay,    a   =    1,    c   =    0.1 
Input    frequency   from    3.5hz    to    3.4   hz. 
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FIGURE  C-15:   1/11  to  1/9  subharmonic  oscillation.   Input 

frequency  from  4.3hz  to  4.2hz.   a  =  1,  c  =  0.1 
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FIGURE  C-17 


Signal  waveform  of  system  with  two  nonlinearities 
C  =  0.387. 

(a)  Input  =  9Sin  3.44TTt  (d)  Output  of  N,  f(t) 

(b)  Error  Signal,  E      (e)  input  to  N?,  e. 

(c)  Input  to  Nx,  e       (f)  System  Output  1 
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FIGURE  C-18:   Singal  waveforms  for  1/5  subharmonic  oscil- 
lation.  Input  =  95in44  t;  K  =50,  k  =0.04. 
r  m    '   t 

(a)  Input  signal  (d)  Input  to  second 

(b)  Error  signal  nonlinearity 

(c)  Input  to  first  nonlinearity   (e)  System  output 
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FIGURE  C-22 


1/3  subharmonic  oscillation  fading  to  funda- 
mental oscillation  due  to  change  of  input 
frequency  from  1.3  to  1.4  hz.   Saturation 
and  backlash:  K   =  0.13,  K  =  100. 
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FIGURE    C-23 


Signal  waveform  showing  change  of  subharmonic 
mode  from  1/4  to  1/5. 


(a)  Input  =  21Sin3.27rt   (d) 

(b)  Error  signal        (e) 

(c)  Input  to  N, 


Input  to  N~ 
System  Output 


101 


LIST  OF  REFERENCES 


1.  Douce,  J.  L.,  Discussion  of  a  Paper  by  Ogata,  Trans. 

ASME,  vol.  80,  no.  8,  p.  1808,  November  1958. 

2.  Gelb,  A.  and  Vander  Velde,  W.  E. ,  Multiple- Input 

Describing  Function  and  Nonlinear  System  Design, 
McGraw-Hill  Book  Company,  1968. 

3.  Gibson,  J.  E.  and  Sirdhar,  R.,  A  New  Dual- Input 

Describing  Function  and  an  Application  to  the  Sta- 
bility of  Forced  Nonlinear  Systems,  Trans.  AIEE, 
pt.  II,  Appl.  Ind.,  p.  65-70',  May  1963. 

4.  Han,  K.  W. ,  Analysis  and  Design  of  Nonlinear  Systems 

by  Critical  Curve  and  Critical  Surf ace  Methods , 
Project  Report  in  connection  with  Nonlinear  Control 
Systems  Course,  Naval  Postgraduate  School,  1962. 

5.  Hayashi,  C.,  Nonlinear  Oscillations  in  Physical 

Systems ,  McGraw-Hill  Book  Company,  1964. 

6.  Ogata,  K.,  Subharmonic  Oscillations  of  Nonlinear  Feed- 

back Control  Systems,  Trans.  ASME,  vol.  80,  no.  8 , 
p.  1802-1808,  November  1958. 

7.  Oldenburger,  R.  and  Nichols,  R.,  Stability  of  Sub- 

harmonic  Oscillations  in  Nonlinear  Systems,  Proc.  JACC, 
Minneapolis,  Minn.,  p.  6  75-680 ,  June  1963 . 

8.  Sawaragi,  Y.  and  Akashi ,  H.,  Condition  for  the  Sub- 

harmonic  Oscillation  in  Relay  Control  System's"!  The 
Science  and  Engineering  Review  of  Doshisha  Univer- 
sity,  vol.  3,  no.  1,  p.  19-35,  May  1962. 

9.  Thaler,  G.  J.  and  Pastel,  M.  P.,  Analysis  and  Design 

of  Nonlinear  Feedback  Control  Systems,  McGraw-Hill 
Book  Company,  196  2. 

10.  West,  J.  C.  and  Douce,  J.  L.,  The  Mechanism  of  Sub- 

harmonic  Generation  in  a  Feedback  System,  Proc.  IEEE , 
vol.  102,  part  B,  no.  5,  p.  569-574',  September  1955. 

11.  West,  J.  L.,  Analytical  Techniques  for  Nonlinear  Con- 

trol Systems,  D.  Van  Nostrand  Company,  Inc.  ,  1960 . 

12.  West,  J.  C,  Douce,  J.  L.,  and  Livesly,  R.  K.,  The 

Dual-Input  Describing  Function  and  its  Use  in  the 
Analysis  of  Nonlinear  Feedback  Systems,  J~.    IEEE , 
London,  vol.  B-103,  p.  463-472,  July  1956. 


102 


INITIAL  DISTRIBUTION  LIST 


No.  Gopies 


1.  Defense  Documentation  Center  2 
Cameron  Station 

Alexandria,  Virginia   22314 

2.  Library,  Code  0212  2 
Naval  Postgraduate  School 

Monterey,  California   93940 

3.  Flag  Officer  In  Command,  Philippine  Navy    1 
Roxas  Boulevard,  Manila 

Philippines 

4.  Professor  G.  J.  Thaler  (thesis  advisor)     2 
Department  of  Electrical  Engineering 

Naval  Postgraduate  School 
Monterey,  California   93940 

5.  Department  of  Physical  Sciences  1 
Philippine  Military  Academy 

Fort  Del  Pilar,  Baguio  City 
Philippines 

6.  Commanding  Officer  1 
Ship  Repair  Yard 

Cavite  Naval  Base,  Cavite  City 
Philippines 

7.  Lt.  (jg)  Hasan  Kocaoglu  1 
Dervis  Ali  mah.  Kefevi  sok.  no.  12 
Karagumruk ,  Istanbul 

Turkey 

8.  Lcdr.  Emerson  C.  Tangan  2 
2671  Int.  11  Juan  Luna 

Gagalangin,  Manila 
Philippines 


103 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  -R&D 

(Security   c  las  silic  ation  ot  title,    body  ot  abstract  and  indexing  annotation  must  be  entered  when   the  overall  report  Is   classified) 


omoiNA  TING    activity   (Corporate  author) 

laval    Postgraduate    School 
Monterey,    California      93940 


2*.  REPORT    SECURITY    CLASSIFICATION 


UNCLASSIFIED 


2fc.     GROUP 


REPOR  T     TITLE 


Investigation   of   Subharmonic   Oscillations    in  Nonlinear   Second   Order 
lystems  . 


DESCRIPTIVE   NOTES  (Type  ot  report  and.  inclusive  dates) 

Master's    Thesis,    December    19  70 


AU  THORISI  (First  name,   middle  initial,   last  name) 

Emerson   Carag   Tangan 


REPOR  T    D A  TE 


December    1970 


7a.     TOTAL    NO.    OF    PAGES 


105 


7b.    NO.    OF    RE  FS 


12. 


a.    CONTRACT    OR    GRANT    NO 


6.    PROJEC  T    NO 


9a.    ORIGINATOR'S    REPORT    NUMBER(S) 


9b.   OTHIIR    REPORT   NO(S)  (Any  other  numbers    that  may  be  assigned 
this  report) 


0      DISTRIBUTION    STATEMENT 


'his  document  has  been  approved  for  public  release  and  sale;  its 
distribution  is  unlimited. 


1       SUPPLEMENTARY    NOTES 


12.    SPONSORING   MILI  T  ARY    ACTIVITY 


Naval  Postgraduate  School 
Monterey,  California   93940 


3.  ABSTRACT 


A  nonlinear  control  system  sometimes  oscillate  at  a  frequency 
that  is  an  integral  submultiple  of  the  driving  frequency  in  re- 
sponse to  a  sinusoidal  input.   Such  a  response  is  undesirable. 
An  analog  computer  study  is  undertaken  to  investigate  the  occur- 
rence of  subharmonic  oscillation  in  second  order  systems  with 
nonlinearities  characterized  by  saturation  and  hysteresis  effects. 


dd  ,fr:..i473 

5/N    0101 -807-681 1 


(PAGE     1 ) 


104 


UNCLASSIFIED 


Security  Classification 


A-SK08 


UNCLASSIFIED 


Security  Classification 


key    wo  R  OS 


Dual- Input  Describing  Function 
Subharmonic  Resonance 
Limit  Cycle 


,Fr:..1473  <back> 

S/N    0101  -807-6821 


UNCLASSIFIED 


105 


Security  Classification 


a- 3  I  409 


Syroeuse,  N.  Y, 
Stockton,  Colif.  ; 


rangan  ••—77 

subharmonic  oscilla- 
st'°nS  '"  "°n'.-nea 
S6COnd  °r<^  systems. 


Thesis 
T136 


1  9 

—  -  ,v 

Tangan 

Investigation  of 
subharmonic  oscilla- 
tions in  nonlinear 
second  order  systems. 


61 J  7 


thesT136 

Investigation  of  subharmonic  oscillation 


3  2768  001  01379  0 

DUDLEY  KNOX  LIBRARY 


